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Abstract. Two completely new algorithms for generating permutations, shift-cursor algorithm 
and level algorithm, and their efficient implementations are presented in this paper. One imple- 
mentation of the shift cursor algorithm gives an optimal solution of the permutation generation 
problem, and one implementation of the level algorithm can be used to generate random per- 
mutations. 
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1. Introduction 

Permutation generation has a history which dates back to the 1650s. After the elec- 
tronic computer was invented dozens of algorithms have been published which generates 
all the permutations of S n by a computer. In 1977 Robert Sedgewick [3] wrote a survey 
of this problem and compared nearly all permutation generation methods at that time. 
As a basic problem of computer science, it is still an interesting problem to be studied. 

Sedgewick [3] concluded in his survey that his implementation of Heap's method 
was the fastest permutation generation algorithm for most computers, and it can be 
coded that most permutations are generated using only two store instructions. This 
is a lower bound of permutation algorithms. Of course it can not be exceeded but 
algorithms which run as fast as it may be written. 

Two completely new algorithms for generating permutations, shift- cursor algorithm 
and level algorithm, are presented in this paper and we give some implementations of 
them. Like Sedgewick's implementation of Heap's method, one of the implementations 
of our first algorithm comes infinitely close to the theoretically optimal solution of the 
problem, that is, 2n! store instructions. Also one implementation of the other algorithm 
can be used to generate random permutations. 

2. The Algorithms 

Our two algorithms are strongly related. Specifically, each number in a "level" permu- 
tation gives the level of the shift cursor in the corresponding permutation. 

First, we discuss the shift cursor algorithm which generates the permutations in S n 
beginning with the identity permutation 123 • • • n. The algorithm uses the first number 
1 in 123. ..n to partition the n! permutation in n blocks of (n — 1)! permutations as 
follows. The algorithm first generates (n— 1)1 permutations with cursor 1 in I s * position, 
the next (n — 1)! permutations with cursor 1 in the 2 n< ^ position, and so on, until the 
cursor 1 is in the n position. In other words, the cursor partitions the n! permutations 
into n blocks, and in each block the cursor is in a fixed position. We define the level 
of the cursor to be the number of blocks determine by the cursor. Thus the level of 
cursor 1 is n. 



The algorithm, while generating the permutations with the cursor 1 in a fixed 
position, uses a new cursor to partition each block of (n — 1)! permutations into n — 1 
smaller blocks and in each smaller block the new cursor is in a fixed position. To find 
the cursor of a smaller block, we first ignore the cursors of the larger blocks from the 
permutations, and then use the first number of the first permutation in the block as the 
cursor of the smaller block. For example, all the level n — \ cursors are the first numbers 
other than cursor 1 in the first permutation of each smaller block. This method can be 
used until the new blocks contain only one permutation and then all n\ permutations 
are generated. 

Consider, for example, the case n = 4 whose permutation sequence appears in Figure 
1. The sequence begins with the permutation 1234. The algorithm first uses 1 as the 
cursor of all 4! permutations, and we get the following 4 blocks (where * stands for an 
unknown number): 



A new cursor is then used in each block to partition the 3! = 6 permutations into 3 
smaller blocks. For example, in the first larger block the first permutation is 1234, so 
we use 2 as the cursor of the first block to get the following 3 smaller blocks: 



12** 




1*2* 




1**2 


12** 




1*2* 




1**2 



In the first block above, where the first permutation is 1234, the first number, other 
than 1 and 2, is 3, so we use 3 as the cursor of the block. Then we get 2 smaller blocks: 
123* and 12*3. Both of the blocks contains only one permutation and we put the last 
number 4 in the permutations to get 1234 and 1243. We continue in this way to obtain 
the other permutation in the first block ending with the permutation 1432. 

In the second block, the first permutation will be 4132. This comes from the previous 
permutation 1432 (the last permutation in the first block) by interchanging the cursor 
1 with the next number in the permutation. Since 4132 is the first permutation in the 
second block, we use the first number 4 as the cursor of the second block to obtain the 
following 3 smaller blocks: 



41** 




*14* 




*1*4 


41** 




*14* 




*1*4 



Continuing as above, we get all the permutations in the second block. We then get 
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all n! = 4! = 24 permutations which appear in Figure 1. 

Observe that, in each permutation, each number in the permutation is a cursor of 
a certain block, and in a given permutation the levels of the cursors are different from 
each other. Thus the cursor levels form a permutation. Our level sequence is obtained 
from the sequence of cursor level permutations by interchanging 1 and n, 2 and n — 1, 
and so on. 

Figure 1 shows, for n = 4, the shift cursor sequence, the levels of the cursors, and 
the level sequence. We also include the well-known Johnson- Trotter sequence [2] [4] for 
reference, (n = 4). See Figure 1. 
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Figure 1. 
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3. Implementation of the Level Algorithm - Also can be 
Used to Generate Random Permutations 

If N is so large that we could never hope to generate all permutations of N elements, it is 
of interest to study methods for generating "random" permutations of N elements. The 
following algorithm implements level algorithm by using a number between 1 to N\ to 
generate a permutation of level algorithm. Since we can generate a "random" number 
between 1 to N\ first and use the following algorithm to generate a permutation of level 
algorithm, the following algorithm can also be used to generate "random" permutations. 

Suppose the function divide(int a, int b) executes a/b and put the quotient to b, 
and the remainder to a and 0! = 1. We can write implementation of level sequence 
using the following algorithm: 

for i=l to n! step 1 do 

begin 

for j=l to n step 1 do 
begin 

int m=(n-j) ! 

divide (i, m) 

if i/0 then 

put j to the (m+l) tl:i empty position 
else 

put j to the m empty position 
end 

output one permutation 
end 

For example, suppose n = 4, to generate the 15*^ permutation we first calculate 
15/3! = 15/6 = 2 • • • 3, the quotient is 2 the remainder is 3, then we put number 1 in 
the 3 rci position, we get Secondly we calculate 3/2! = 3/2 = 1 • • • 1, the quotient 

is 1 the remainder is 1, then we put number 2 in the 2 nc ^ empty position, we get *21*. 
Third we calculate 1/1 ! = 1 • • • 0, the quotient is 1 the remainder is 0, then the number 
3 is in the 1 st empty position, we get 321*. Finally we calculate 0/0!=0/l=0...0, the 
quotient is the remainder is 0, then the number 4 is in the 1 st empty position, we get 
the permutation 3214 at last. 

4. Implementation of Shift Cursor Algorithm - an Optimal 
Solution of Permutation Generation Problem 

It is can be noticed that shift cursor sequence can generate next permutation by trans- 
position. And the sequence has the following properties: 

Property 1. A level j cursor remains the same position in continuous (j — 1)! 
permutations. 

Property 2. // ignore the cursors which level > j from a sequence of n! permu- 
tations then the sequence turns to n\/j\ blocks with length j, and the rule of change of 
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each blocks is the same. 

We call Property 1 as the local stability of the sequence, and call Property 2 as the 
rule stability of the sequence. 

Property 1 ensures that when generating permutations using shift cursor algorithm, 
cursors with higher levels can be placed in fixed positions for a certain part of the 
sequence. 

Property 2 ensures that for the certain part of the sequence, if higher level cursors 
have already been placed in their positions, a predetermined rule can be used to place 
other cursors in their right positions. 

Clearly [1], there are two types of operations needed for generating a new permu- 
tation from the previous permutation: 

1. Operations needed for actually executing changes from the previous permutation. 

2. Operations needed for reaching decisions about which type 1 operations are re- 
quired. 

In theory, any permutation generating algorithm should contain the two types of 
operations and only the two types of operations, and a theoretically optimal algorithm 
for the problem should has the following factors: 

1. All type 1 operations should be change only two numbers of the permutation. 

2. The number of type 2 operations should be limited to reach zero, or much smaller 
than type 1 operations. 

By using the two properties of shift cursor sequence, algorithms which have the above 
two factors can be written. 

Suppose we have already written an algorithm for generating permutations of shift 
cursor algorithm. For a certain part of the sequence, such as a sequence of continuous 
j\ permutations, we can place the cursors which level > j in their places first, and 
these cursors remain the fixed positions in the j! permutations (property 1). Then by 
using a given rule, we can place other cursors directly (property 2), without any type 
2 operations. 

For example, suppose we are using the algorithm to generate 5! permutations 
(n = 5), and suppose j = 3. Starting from 12345, the level of number 1 is 5, the level 
of number 2 is 4. Then we fix number 1 in the 1 st position, fix number 2 in the 2 nc ^ 
position. Afterwards, by using a given rule, we can place the numbers 3, 4, 5 directly 
in their positions to generate permutations. The rule to place the three numbers is 
the same as the rule to generate 3! permutations (n = 3) using shift cursor algorithm. 
When we using the algorithm to generate 3! permutations (n = 3), the sequence is as the 
following: 123, 132, 312, 213, 231, 321. In this case we name 3 the No.l element, name 
4 the No. 2 element, name 5 the No. 3 element. In the above sequence, No. 2 element and 
No. 3 element be changed firstly, so we change 4, 5 first, then the next permutation to 
12345 is 12354. And No.l element and No. 3 element changed afterwards, so we change 
3, 5 secondly, then the next permutation to 12354 is 12534, and so on. And finally we 
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get the 6 permutation 12543. Some type 2 operations is needed now to determine 
which cursor's level > j and where should we place these cursors. After that, we can 
generate next 3! — 1 = 5 permutations by directly change two numbers. This algorithm 
can be described as following: 
Algorithm 1: 

1. Starting from 1, 2, 3, ... ,n. 

2. If the level order of the cursors is 1, 2, 3, ... , n, stop. 
3. Determine cursors which level>j and their positions. 
4 . Output . 

5. Generate and output next j!-l permutations directly. 

6. Go to step 2. 

When generating permutations using algorithm 1, there are (j\ — 1) / j! permutations 
are generated by directly change two numbers without any type 2 operations. And other 
permutations are generated using additional type 2 operations. 

Suppose when using algorithm 1, the average execution time of type 1 operation is 
t and the average execution time of type 2 operation is T, then the overall execution 
time of algorithm 1 (Tl) is n! *t + n\/j\ *T. Since every type 1 operation of algorithm 1 
is change two numbers, the overall execution time of theoretically optimal algorithm of 
permutation generation (TO) is n! * t, and T1/T0= 1 + Tj (j\ * t) > 1. So with j grows 
larger T1/T0 could infinitely close to 1, which means algorithm 1 could infinitely close 
to the theoretically optimal solution of the problem. 

Robert Sedgewick's implementation of Heap's method is very similar to algorithm 1, 
the only difference between them is the termination condition (step 2) of the algorithms, 
and they both have the property 1 and property 2. And it is the two properties make 
the two algorithms both can infinitely close to the theoretically optimal solution of the 
problem. 

5. Summary and Conclusions 

Although there are few applications need to generate all n! permutations, permutation 
generation is still an interesting problem to be studied. Further studies are needed on 
this problem, and these studies may change our understanding of permutations in the 
future. 

Comparing with permutation generation, random permutation generation has more 
practical meanings. The implementation of level algorithm gives a new method to 
generate random permutation. And it can be used to generate the permutation which 
in a given position of the level sequence. 

Because of the two important properties, the shift cursor algorithm has an efficient 
implementation, and in theory both the implementation and Sedgewick's implemen- 
tation of Heap's method could infinitely close to the theoretically optimal solution of 
the problem, which is 2n! store instructions. Comparing with Heap's algorithm and 
another famous algorithm - Johnson- Trotter algorithm, the shift cursor algorithm is 
easier to understand than Heap's algorithm and has a more efficient implementation 
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than Johnson- Trotter algorithm. And it's more suitable to be written in textbooks. 
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